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Abstract 

Thermo dynamical properties of canonical and grand-canonical ensembles of the half-filled 
two-site Hubbard model have been discussed within the framework of the nonextensive 
statistics (NES). For relating the physical temperature T to the Lagrange multiplier /3, 
two methods have been adopted: T = l/ksfl in the method A [Tsallis et al. Physica A 
261 (1998) 534], and T = c q /k B (3 in the method B [Abe et al Phys. Lett. A 281 (2001) 
126], where k B denotes the Boltzman constant, c q = J2iPi, Pi t ne probability distribution 
of the ith state, and q the entropic index. Temperature dependences of specific heat and 
magnetic susceptibility have been calculated for 1 = q = 2, the conventional Boltzman- 
Gibbs statistics being recovered in the limit of q = 1. The Curie constant T q of the 
susceptibility in the atomic and low-temperature limits (t/U — > 0, T/U — > 0) is shown 
to be given by T g = 2 q 2 2( - q ~ 1 ' > in the method A, and T g = 2 q in the method B, where t 
stands for electron hoppings and U intra-atomic interaction in the Hubbard model. These 
expressions for T q are shown to agree with the results of a free spin model which has been 
studied also by the NES with the methods A and B. A comparison has been made between 
the results for canonical and grand-canonical ensembles of the model. 
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1 INTRODUCTION 



In the last several years, there is an increased interest in the nonextensive statistics (NES), 
which was initially proposed by Tsallis [1, 2, 3]. This is because the standard method 
based on the Boltzman-Gibbs statistics (BGS) cannot properly deal with nonextensive 
systems where the physical quantity associated with N particles is not proportional to N 
[4]. The nonextensivity has been realized in three different classes of systems: (a) systems 
with long-range interactions, (b) small-scale systems with fluctuations of temperatures or 
energy dissipations, and (c) multifractal systems [3, 5]. In a gravitating system with the 
long-range interaction, which is a typical case (a), the specific heat becomes negative [6]. 
A cluster of 147 sodium atoms, which belongs to the case (b), has been reported to show 
the negative specific heat [7]: note that the specific heat is never negative in the canonical 
BGS [8]. 

Tsallis [1] [2] has proposed a generalized entropy in the NES defined by 

s q = k B , (1) 

where ks denotes the Boltzman constant, pi the probability distribution of the system in 
the ith configuration, and q the entropic index. The entropy of BGS, Sbg, is obtained 
from Eq. (1) in the limit of q — 1, as given by 

Si = Sbg = -k B ^Pi hi Vi- (2) 

i 

The nonextensivity in the Tsallis entropy is satisfied as follows. Suppose that the total 
system is divided into two independent subsystems with the probability distributions, 
p^ and pf\ The total system is described by the factorized probability distribution 
Pij = p^ pf \ The entropy for the total system given by Eq. (1) satisfy the relation: 

S q = SU+Sp + (^)sV>Sf\ (3) 

where S^ stands for the entropy of the kth subsystem. Equation (3) shows that the 
entropy is extensive for q — 1 and nonextensive for q ^ 1: the quantity | q — 1 | expresses 
the measure of the nonextensivity. The new formalism has been successfully applied to a 
wide range of nonextensive systems including physics, chemistry, mathematics, astronomy, 
geophysics, biology, medicine, economics, engineering, linguistics, and others [5]. 
The current NES, however, is not complete, having following unsolved issues. 
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(i) A full and general understanding of the relation between the entropic index q and 
the underlying microscopic dynamics is lacking. The index q is usually obtained in a 
phenomenological way by a fitting of experimental or computational available data. It 
has been reported that the observed velocity distribution of galaxy clusters significantly 
deviates from BGS distribution, which may be fitted well by the NES distribution with 
q ~ 0.23 [9]. The index q in self-gravitating systems has been shown to be q = 0.60 ~ 0.92 
for n = 3 ~ 20 with the use of the relation: n — 1/(1 — q) + 1/2 for the entropic index 
q and the polytrope index n given by P(r) oc p(r) 1+1 / n , where P{r) is pressure and p{r) 
mass density [10]. However, the index q has been successfully determined a priori in 
some cases [3]. In one-dimensional dissipative map, the index is given by q = 0.2445 [11] 
from its scaling property of dynamic attractors. In nanometric systems consisting of N 
noninteracting particles, the NES distribution is shown to arise from fluctuating (3 whose 
distribution is given by the x 2 ( or H distribution, leading to q — 1 + 2/N [12, 13, 14]. 

(ii) The second issue is that it is not clear how to relate the physical temperature T to 
the introduced Lagrange multiplier (5. So far two methods have been proposed: 

T = jL, (method A) (4) 

= (method B) (5) 

where c q = J2iPi- The method A proposed in Ref. [2] is the same as the extensive BGS. 
The method B is introduced so as to satisfy the zeroth law of thermodynamical principles 
and the generalized Legendre transformations [15]. It has been demonstrated that the 
negative specific heat of a classical gas model which is realized in the method A [16], is 
remedied in the method B [15]. Results calculated with the use of the two methods have 
been compared in self-gravitating systems [17]. 

In this paper, we wish to apply the NES to the Hubbard model, which is one of the 
most important models in condensed-matter physics (for a recent review on the Hubbard 
model, see Ref. [18]). The Hubbard model consists of the tight-binding term expressing 
electron hoppings and the short-range interaction between two electrons with opposite 
spins. The model provides us with good qualitative description for many interesting 
phenomena such as magnetism, electron correlation, and superconductivity. In particular, 
the model has been widely employed for a study of metallic magnetism. In the atomic 
limit where the electron interaction is much larger than electron hoppings, the Hubbard 
model with the half-filled band reduces to a local-spin model such as the Heisenberg and 
Ising models. Despite the simplicity of the Hubbard, however, it is very difficult to obtain 
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its exact solution. In order to get a reasonable solution, various technical methods have 
been developed [18]. It is possible to employ a small number of lattice models to obtain an 
analytical solution. The two-site Hubbard model is employed as a simple model which can 
be analytically solved. Thermodynamical and magnetic properties of the two-site model 
have been studied within the BGS [19, 20, 21]. Although the two-site Hubbard model is a 
toy model, it is an exactly solvable quantum system. Actually this model may be adopted 
to describe interesting phenomena in real systems such as organic charge-transfer salts 
with dimerized structures [19]-[22]. Because the two-site Hubbard model is considered to 
belong to the case (b) of small-scale systems mentioned above, it is worthwhile to apply 
the NES to the model in order to investigate nonextensive effects on its thermodynamical 
and magnetic properties, which is the purpose of the present paper. 

The paper is organized as follows. After discussing the NES for grand-canonical en- 
sembles (GCE) of the two-site Hubbard model in Sec. 2, we will calculate the specific 
heat, susceptibility and the Curie constant of the susceptibility in the atomic limit calcu- 
lated by the NES with the two methods A and B for the T — (3 relations given by Eqs. 
(4) and (5). In Sec. 3, we present calculations for canonical ensembles (CE) of the model, 
which are compared to those for GCE. The final Sec. 4 is devoted to discussions and 
conclusions. In appendix A, free spin systems have been discussed also by using the NES 
with the two methods A and B. 

2 Grand-canonical ensembles 

2.1 Entropy, energy and free energy 

We consider GCE of the half-filled, two-site Hubbard model given by 



where nj a = aj a aj a , aj a denotes the annihilation operator of an electron with spin a on 
a site j (=1, 2), t the hopping integral, U the intraatomic interaction and h an applied 
magnetic field in an appropriate unit. Eigen values of the system are given by [21] 



2 2 



H = -t^2{a\ a a 2a + a\ a a la ) + U ^ n n n n - hJ2(n n - n n ), 



(6) 




= ±t±h, 



0, ±2h, U, U/2 ± A, 



for i — 1 (rii = 0) 
for % = 2 to 5 (rii = 1) 
for i = 6 to 11 (rii = 2) 



(7) 



(9) 
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U±t±h, 
2U, 



for % — 12 to 15 (rii = 3) 
for i = 16 (ni = 4) 



(10) 
(11) 



where A = y f/ 2 /4 + At 2 , and expresses the number of electrons in the ith state. The 
grand-partition function in BGS, S_bg, is given by {ks = 1 hereafter) [21] 



= 1 + 4 cosh(/i/T) cosh(t/T) e^ /T 

+ [1 + 2 cosh(2/i/T) + e"^ + 2 e~ u/2T cosh(A/T)]e 2 ^ T 
+4 cosh(/i/T) cosh(t/T) e-( u - 3 ^ T + e -2(^- 2 ")/ r . 



'12) 



(13) 



Here Tr stands for trace, N = Y%=i J2 a n icy> an d the chemical potential is p = U/2 inde- 
pendent of the temperature in the half-filled case where the number of total electrons is 
N e = 2. By using the standard method in the BGS, we can obtain various thermodynam- 
ical quantities of the system [21]. 

Now we adopt the NES in which the entropy S q is given by [1] [2] 

_ (Tr(pl)-1\ 



l-q 



(14) 



Here p q denotes the generalized density matrix, whose explicit form will be determined 
shortly [Eq. (18)]. We impose the three constraints given by 



Tr (Pg) 
Tr (g H) 

Tr [p\) 
Tr (pi N) 



= < H >q= E q , 

= <N>g=N q , 



(15) 
(16) 

(17) 



Tr ($) 

where the normalized formalism is adopted [2]. The variational condition for the entropy 
with the three constraints given by Eqs. (15)-(17) yields 



with 



X„ 



Pq = !T eX Vq 
A q 



Tr exp 



£-) (H — pN — E q + pN q ) 



- (H-pN-E q + pN q ) 



(18) 



c g = Tr(M)=X, 



(19) 
(20) 
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where exp Jx) = [1 + (1 — q)x]~ is the generalized exponential function. Lagrange 
multipliers f3 and p relevant to the constraints given by Eqs. (16) and (17) are given by 
the relations: 

* = §■ (2i) 

n - _I^L (22) 
The entropy 5 9 in Eq. (14) is expressed by 

S "={ X [ q -~ 1 )= ] ^ < 23 ) 

where ln^rr) = (a; 1-9 — 1)/(1 — 9) is the generalized logarithmic function. 

In relating the physical temperature T to the Lagrange multiplier (3, we have adopted 
the two methods [2] [15] : 

T = i (method A) (24) 

= ^. (method B) (25) 

In the limit of q — 1, we get the results obtained in the BGS: p BG = e~^ H ~^ /E BG , 
E BG = Tr(p BG H), N BG = Tr(p BG N), X BG = e ^ E ^-^ BG ) Sbg? Ebg = Tr ( e -P(H-»N)^ 

and S BG = -Tr (p BG \np BG ). 

It is necessary to point out that E q and X q have to be determined self-consistently 
by Eqs. (16)-(20) with p, determined by Eq. (22) for N q = N e and a given temperature 
T because they are mutually dependent. In the half-filled case, however, calculations 
become easier because p = U/2 independent of the temperature. The calculation of 
thermodynamical quantities in the NES generally becomes more difficult than that in 
BGS. In our numerical calculations to be reported in this paper, simultaneous equations 
for E q and X q given by Eqs. (16)-(20) are solved by using the Newton-Raphson method. 
The iteration have started with initial values of E q and X q obtained from the BGS (q = 1). 
Numerical calculations have been made for 1 = q = 2 which is appropriate for nanoscale 
systems[12, 13, 14]. 

Figures l(a)-l(f) show the temperature dependence of the energy E q of GCE calculated 
for h — 0. Bold solid curves in Fig. 1(a), 1(b) and 1(c) show E\ in the BGS calculated 
for U/t — 0, 5 and 10, respectively. The ground-state energy at T = is E x jt = -2.0, 
-0.70156 and -0.38516 for U/t — 0, 5 and 10, respectively. With increasing q value above 
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unity, the gradient of E q is much decreased in the method A, as shown in Figs. l(a)-l(c). 
This trend is, however, much reduced in the method B, as Figs. 1 (d)-l(f ) show. This 
behavior is more clearly seen in the temperature dependence of the specific heat C q , as 
will be discussed shortly [Figs. 3(a)-3(f)]. 

Temperature dependences of the entropy for h — are plotted in Figs. 2(a)-2(f). 
Figures 2(a), 2(b) and 2(c) express S q for U/t — 0, 5 and 10, respectively, calculated by 
the method A, and Figs. 2(d)-2(f) those calculated by the method B. Bold curves denote 
the results for the BGS, where the entropy is quickly increased at low temperature when 
the interaction is increased. When the q value is more increased above unity, S q is more 
rapidly increased at very low temperatures and its saturation value at higher temperatures 
becomes smaller. This behavior is commonly realized in the results calculated by the 
methods A and B. A difference between the two results is ostensibly small because S q 
shows a saturation at low temperatures. 

2.2 Specific heat 

First we consider the specific heat, which is given by 



Because E q and X q are determined by Eqs. (16)-(20), we get simultaneous equations for 
dE q /d[3 and dX q /d{3, given by 




(26) 




(27) 



(28) 



with 



a u = q(3X q q 2 J2 w i q 



(29) 



a 12 = -X- l E q -Pq{q-l)X^Y. 




e-i - jirii - E q + (iN e ), 



(30) 



a 22 = 



G21 = PXl 



2q-l 



€i(€i — E q ), 



(31) 
(32) 
(33) 



Wi = < i | exp g — I — J (H — jiN — E q + /iN q ) \ i >, 
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!-(!-?) 



1-9 



The specific heat is then given by 



(34) 
(35) 



V dry VI - an - ai 2 a 2 i 



(36) 



with 



9/5 
3T 



dX q 
dp 



~P\ 



P 2 



X^ q -P(l-q)X^(dX q /dp), 
bi 



(method A) 
(method B) 



(1 - an - a l2 a 2i y 
(1 - an - ai 2 a 2 i)' 



(37) 
(38) 

(39) 
(40) 



It is worthwhile to examine the limit of q = 1 of Eqs.(29)-(35), which reduce to 

an = pE 1 , 

aw = —X 1 1 E±, 

a 2 i = 0X U 

bi = - < e\ >i + < €i >\ . 



yielding 



with 



d = P\< ef > 1 -<e l >i), 



<Qi>! = Xf 1 e -^-^- El+ ^ Q h 

i 

E x = <ti >!, 



xi = X 



D -f3(£i-^n i -E 1 +fMN e ) 



(41) 
(42) 
(43) 
(44) 

(45) 

(46) 

(47) 
(48) 



The expression for C\ agrees with Cbg obtained in the BGS. 

Figures 3(a)-3(f) show the specific heat of GCE calculated for h — 0. C\ in BGS 
for U/t = shown by the bold solid curve in Fig. 3(a), has a peak at T/t ~ 0.65. 
Figure 3(c) shows that for U/t = 10, this peak splits into two. A lower peak arises 
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from low-lying collective spin-wave-like excitations while higher one from single-particle 
excitations [20] [21]. For intermediate U/t = 5 these two peaks overlap [Fig. 3(b)]. The 
temperature dependences of the specific heat C q calculated with the use of the method A 
for U/t — 0, 5 and 10 are plotted in Figs. 3(a), 3(b) and 3(c), respectively. We note that 
when q is larger than unity, peaks become broader. Figures 3(d), 3(e) and 3(f) show the 
temperatures dependence of the specific heat C q calculated by the method B for U/t — 0, 
5 and 10, respectively. Although general property of the q dependence of the specific heat 
of the method B is similar to that of the method A, the effect of the nonextensivity in 
the method B becomes smaller than that in the method A. 

2.3 Susceptibility 

In the NES, the magnetization induced by an applied field h is obtained by 

m q = < Hi > q , (49) 
leading to the susceptibility given by 

Xq = ~dh Ih=0 ' ^ ' 

= -Ef^(3-'X^Xf\ (51) 

where jXi = — dei/dh and = d 2 E q /dh 2 \h=o et al. With calculations using Eqs. 
(16)-(20), we get simultaneous equations for E^ and X^\ given by 

= anE^ + anX^ + n, (52) 

X® = a 21 E^ + a 22 X^+f 2 , (53) 

with 

h = -2/??i:- 2 E^VL (54) 

h = PUX^^w^nl (55) 

i 

where = 1,2) are given by Eqs. (29)-(31). The susceptibilities Xq i n t ne methods 

A and B are expressed by 

= / +r . x -. (1 _ ail) x = K 
\ 1 - an - 012^21 / a 21 

with the T — (3 relations given by Eqs. (24) and (25), respectively. 
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In the limit of q = 1, Eqs. (54) and (55) reduce to 

fi = -2p(i4) lt (57) 
h = f3 2 X, (/i^, (58) 

leading to the susceptibility in the BGS: 

Xbg = Xi = P<tf >i> ( 59 ) 

where the bracke < • >i is given by Eq. (46). 

The BGS susceptibility for U/t = has a peak at T/t ~ 0.65 as Fig. 4(a) shows. 
With increasing U/t, the magnitude of Xbg is enhanced by the interaction, and its peak 
positon becomes lower [20] [21], as Figs. 4(b) and 4(c) show: the horizontal scale of Fig. 
4(c) is enlarged compared with Figs. 4(a) and 4(b). The temperature dependences of the 
susceptibility Xq calculated by the method A for U/t — 0, 5 and 10 are plotted in Figs. 
4(a), 4(b) and 4(c), respectively. We note that as increasing q above unity, the peak in 
Xq becomes broader. Figures 4(d), 4(e) and 4(f) show the temperature dependence of the 
susceptibility Xq calculated by the method B for U/t — 0, 5 and 10, respectively. Again 
the effect of the nonextensivity in the method B becomes smaller than that in the method 
A. 



2.4 Curie constant in the atomic limit 

The half-filled Hubbard model with t/U C 1 reduces to a local-spin model with the 
superexchange interaction J ~ t 2 /U. In the limit of t/U ^ (atomic limit) for which 
J — > 0, the susceptibility of the two-site Hubbard model in BGS is given by 

2 

Xbg ~ T(l + e^/ 2T )' (6 °^ 
= |, forT«C/ (61) 

= 1. forT>[/ (62) 

Defining the effective, temperature-dependent Curie constant T g (T) by 

T q (T)=T Xq (T), (63) 

we note from Eqs. (61) and (62) that it varies from Ti(0) = 2 for localized moments to 
ri(oo) = 1 for delocalized moments. 
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By using the NES, we get the Curie constant of our model in the low-temperature 
limit (T U), given by 

r,(0) = 2qM q - 1 = 2q2 2{q - 1) , (method A) (64) 

= 2 q, (method B) (65) 

with 

X q = M = 2 2 , (66) 

where M denotes the number of states with the lowest value of e« — /in^. Similarly the 
Curie constant in the high-temperature limit (T ^> U) is given by 

r 9 (oo) = q M q ^ 1 = q A 2{q ~ 1] , (method A) (67) 

= q, (method B) (68) 

with 

X q = Moo = 4 2 , (69) 

where expresses the number of available states in our model. Expressions given by 
Eqs. (64) and (65) are consistent with the results for free spin systems calculated by NES 
with the methods A and B [Eqs. (90) and (91)], whose detail is discussed in appendix A. 

Solid curves in Fig. 5(a) and 5(b) show the temperature dependence of the inversed 
susceptibility l/x g for t = with various q values calculated by the methods A and B, 
respectively. We notice that as increasing the q value, the gradient oil/xq is decreased in 
both the methods. Solid curves in Fig. 5(c) and 5(d) express the temperature dependence 
of Y q calculated by the methods A and B, respectively. We note that r x is 2 at T/U ~ 
and approaches 1 at T/U ~ 10, as Eqs. (61) and (62) show. T 2 (T) in the method A is 16 
at T/U ~ and it is rapidly increased as increasing the temperature with the maximum 
at T/U ~ 4, above which it is decreased. In contrast, T 2 (T) in the method B starts from 
4 at T/U ~ 0, and approaches 2 at T/U = 10 with a peak at T/U ~ 0.2. 

We note that in Figs. 5(a)-5(d) that the Curie constant is increased with increasing 
q. This is more clearly seen in Figs. 5(e) and 5(f), where we plot the q dependences of 
T ? (T) calculated by the methods A and B, respectively. Circles in Fg. 5(e) show T ? (T) at 
T = 0.02 calculated by the method A, which nicely obeys the relation T q (0) = 2q 2 2{q ~ 1) 
[Eq. (64)]. 

On the contrary, r g (oo) in the model A is expressed by squares in Fig. 5(e) where 
the chain curve denotes r 9 (oo) = q A 2 ^^ [Eq. (67)]. A disagreement between the chain 
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curve and the result of T/U = 10 arises from a fact that the temperature of T/U = 10 
does not correspond the high-temperature limit, which is realized at T/U > 100. Circles 
in Figure 5(f) show the q dependence of T q (T) at T/U = 0.02 calculated by the method 
A, which follows the relation T g (0) = 2q given by Eq. (65). On the contrary T g (T) at 
T/U = 10 calculated by the method B is plotted by squares in Fig. 5(f), approximately 
following the chain curve given by r 9 (oo) = q [Eq. (68)]. 



3 Canonical ensembles 

In the previous section, we have discussed thermodynamical properties of GCE of the 
two-site Hubbard model. It is straightforward to extend our study to CE of the model. 
The canonical partition function Z B g in the BGS is given by [19] 



Z BG = l + 2 cosh(2/i/T) + e- u/T + 2 e~ u/2T cosh(A/T). 



(70) 



The generalized density matrix p q for CE in the NES may be determined by the 
variational condition for the entropy given by Eq. (14) with the two constraints given by 
Eqs. (15) and (16) 



with 



X„ 



Tr exp 



(71) 

(72) 
(73) 



where f3 stands for the Lagrange multiplier and E q is given by Eq. (16). The T — (3 
relations in the methods A and B are given by Eqs. (24) and (25), respectively. The 
specific heat is expressed by Eqs. (29)-(40) but with a 12 , Wi and X q given by [23] 



aia = -X; 1 E q -(3q(q-l)Xf 3 Y, w ?~ 1 ^-E q ), 



Wi 

X n 



1 - (1 - q) 



(74) 

(75) 
(76) 



where a sum over i {— 6 to 11) is only for n« = 2 in Eq. (9). The susceptibility is similarly 
expressed by Eqs. (54)-(56) with 012, Wi and X q given above. The Curie constants with 
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t/U = in the low-temperature limit are given by Eqs. (64) and (65) in the methods A 
and B, respcitively, while those in the high-temperature limit are given by Eqs. (67) and 
(68) but with q -> (4/3) q and = C| where CI = n\/k\(n - k)\. 

Temperature dependences of the specific heat C q of CE calculated for h = are plotted 
in Figs. 6(a)-6(f). Bold curves in Figs. 6(a)-6(c) show the temperature dependence of 
C B g in the BGS. For U/t = 0, C B g calculated for CE is similar to that for GCE shown in 
Fig. 3(a). For U/t = 10, however, the result for CE is rather different from that for GCE 
shown in Fig. 3(c): the latter has a broad peak at T/t > 1 while the former not. The 
high-temperature peak in Cbg of GCE is responsible to charge fluctuations whose effect 
is included in the GCE treatment but not in the CE one. The temperature dependences 
of the specific heat in the NES calculated by the method A are plotted in Figs. 6(a)-6(c), 
whereas those calculated by the method B are shown in Figs. 5(d)-5(f). 

Figures 7(a)-7(c) show the temperature dependence of the susceptibility in the NES 
calculated by the method A. Results calculated by the method B are plotted in Figs. 
7(d)-7(f). The difference between the results of CE and GCE is small. This is because 
the susceptibility is not sensitive to charge fluctuations. Effects of the nonextensivity on 
C q and Xq i n the method B become smaller than those in the method A for CE, just as 
for GCE. 

Although results calculated for CE and GCE coincide for infinite systems, it is not 
the case for finite ones. The CE-method is appropriate for an analysis of real nanoscale 
materials where the number of total electrons is conserved. In contrast, the GCE-method 
is expected to be better in guessing properties of infinite systems from results of finite 
systems [20]. This is understood as follows. Suppose that an infinite system is divided 
into many segments. Charge fluctuations exist between adjacent segments, which are 
taken into account in the GCE statistics. 

4 Discussions and Conclusions 

Recent progress in atomic engineering makes it possible to create nanoscale materials 
with the use of various methods (for reviews, see Ref. [24]). Nanomagnetism shows 
interesting properties different from bulk magnetism. Nanoclusters consisting of transi- 
tion metals such as Fe N (N= 15-650) [25], Co^v (N=20-200) [26], and Ni N (N=5-740) 
[27] have been synthesized by laser vaporization and their magnetic properties have been 
measured, where N denotes the number of atoms per cluster. Magnitudes of magnetic 
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moments are increased with reducing N [27]. It is shown that magnetic moments in 
Co monatomic chains constructed on Pt substrates are larger than those in monolayer 
Co and bulk Co [28]. Recently Au nanoparticles with average diameter of 1.9 nm (in- 
cluding 212 atoms), which are protected by polyallyl amine hydrochloride (PAAHC), are 
reported to show ferromagnetism while bulk Au is diamagnetic [29]. This is similar to 
the case of gas-evaporated Pd fine particles with the average diameter of 11.5 nm which 
show the ferromagnetism whereas bulk Pd is paramagnetic [30]. The magnetic property 
of 4 Ni atoms with the tetrahedral structure in magnetic molecules of metallo-organic 
substance [Moi 2 3 o(a*2 - OH)ioH 2 {Ni(H 2 03)}4] • 14H 2 has been studied [31]. Their 
temperature-dependent susceptibility and magnetization process have been analyzed by 
using the Heisenberg model with the antiferromagnetic exchange couplings between Ni 
atoms [31]. Similar analysis has been made for six-, eight-, ten- and twelve-membered Fe 
atoms in a family of magnetic molecules [32]. 

It is well known that the Hubbard model is more suitable for a study on bulk mag- 
netism of transition metals than the Heisenberg or Ising model which is best applied to 
rare-earth-metal magnets like Gd [18]. This is true also for transition-metal nanomag- 
netism. It is interesting to make a comparison of our result obtained for the two-site 
Hubbard model with experimental data of nanoclusters consisting of two transition-metal 
atoms. Our NES calculations have shown that the results for q = 2, which is realized 
for a two-atom molecule if we assume q = 1 + 2/N [12, 13, 14], are rather different from 
those for q — 1 (BGS). Unfortunately ideal clusters including only two transition-metal 
atoms are not available at the moment. Some charge-transfer salts like tetracyanoquin- 
odimethan (TCNQ) with dimerized structures, have been analyzed by using the two-site 
Hubbard model within the BGS [21]. Their susceptibility and specific heat were studied 
by taking into account the effect of the interdimer hopping, which is smaller than the in- 
tradimer hopping, by a perturbation to the two-site Hubbard model [21]. For an analysis 
of nanochains including multiple dimmers within the NES, it is necessary to know the 
size dependence of the q value appropriate for them [33]. 

For a better understanding of nanomagnetism, we are requested to perform NES cal- 
culations of the iV-site Hubbard model with a larger N in one, two and three dimensions. 
The number of eigen values is 4^ for the GCE and C™ for the CE with the half-filled 
electron occupancy. When the orbital degeneracy is taken into account for a more realistic 
description, these numbers become 4 ND and Cffjj' for GCE and CE, respectively, where 
D (= 2L + 1) denotes the orbital degeneracy for the orbital quantum number L (L = 2 
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for 3d transition metals) [34]. Realistic NES calculations for larger nanosystems become 
much laborious. We hope that by extending the current engineering technique, it would 
be possible to synthesize ideal, small nanoclusters consisting of transition-metal atoms. 

To summarize, we have discussed thermodynamical and magnetic properties of GCE 
and CE of the two-site Hubbard model within the frame work of the NES. This is the first 
application of NES to thermodynamical and magnetic properties of the Hubbard model, 
as far as the author is aware of. The temperature dependences of the specific heat and 
susceptibility calculated by the two methods A and B for T — (3 relation, are qualitatively 
the same, and they change significantly when q deviates from unity. The two methods, 
however, yield a quite different q dependence of the Curie constant, which is demonstrated 
for the atomic-limit case, and which is consistent with the result for free spin models as 
discussed in appendix A. Although the result of the method B is considered to be more 
reasonable than that of the method A, we reserve to conclude which of the two methods 
is correct for the two-site Hubbard model at this stage. It is necessary to further develop 
the NES theory in order to clarify the two issues (i) and (ii) raised in the introduction. 
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Appendix A NES for a free spin model 

We consider Hamitonian expressing the N Ising spins with S = 1/2, given by 

H = -ghJ2S ZJ = -ghS z , (77) 

where g{— 2) denotes the g factor and h an applied magnetic field in an appropriate units. 
Eingen states are classified by m s = S z , which runs from —N/2 to N/2. The multiplicity 
of the state specified by m s is given by [35] 

AH 

/(m " N) = W 2 + m.)l(Af/2-m.)| - (?8) 



which satisfies the sum rule: 



N/2 

]T f(m s ,N)=2 N = M. (79) 

m s =-N/2 
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Hereafter we specify the state with the index i, which runs from 1 to M. In the case of 
N = 2, for example, we get rrii = —1 for % = 1, rrii = for i — 2,3 and m 8 = 1 for i = 4. 
By employing the NES for CE as in Sec. 3, we get the probability density given by 



Pi 



X' 



Wi 



cxp„ 



— J (hrrii + E q 



i 

E q = X^Y.hrrii exp ? 



0\ x<? 
1 (hrm + E q ) 



(80) 
(81) 

(82) 
(83) 

(84) 



The thermal average of the magnetization is given by 



m q = g 

from which we get the susceptibility given by 

dm r , , 



- i Pi ' 

E*p1 



Xq 



h=0, 



9 



dh 

2 ( qf3 



1 2q-\ 2 



In deriving Eq. 



we have adopted Wi — 1 and X„ — 2 , and 



(85) 

(86) 
(87) 

(88) 



N/2 



2~ N ]T f(m.,N) 



iV 

T' 



(89) 



m s = -N/2 

for (f3/c q )(h rrii + E q ) <^ 1 with h — > 0. When we adopt the method A given by T = 1//? 
[Eq. (24)], Eq. (88) becomes 



X q 



ZJL 
T ' 



with the Curie constant r„ given by 



NqM q ~ l =Nq2 N(q - 1 \ 



(90) 



(91) 



which agrees with the result of Ref. [36] previously obtained by using the Tsallis' nor- 
malized scheme [2]. Equation (91) has an anomalous exponential dependence on N, and 
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shows dark magnetism: the apparent number of spins is larger than the actual one [36, 37]. 
On the contrary, when we adopt the method B given by T = c q / (3 [Eq. (25)], the Curie 
constant becomes 



Nq, 



(92) 



whose N dependence seems reasonable. 
By using the T — (3 relation given by 



1 (3 



(93) 



T c q + {l-q)(3E q ' 



Reis et al. [35] have obtained the result same as Eq. (91). This is due to a fact that with 
E q = for h = 0, Eq. (92) reduces to T = c q / (3 which is nothing but the T — (3 relation 
in the method B. 

When we adopt the relation: q — 1 + 2/N obtained for nanoscale systems [12, 13, 14], 
Eqs. (90) and (91) become 




(method B) 



(method A) 



(94) 



(95) 



Equation (94) is consistent with the result of the Ising model with A" — > oo. 



18 



References 



[1] C. Tsallis, J. Stat. Phys. 52 (1988) 479. 

[2] C. Tsallis, R. S. Mendes and A. R. Plastino, Physica A 261 (1998) 534. 
[3] C. Tsallis, Physica D 193 (2004) 3. 

[4] There are some systems with long-range interactions, whose energy and free energy 
are not proportional to N while the entropy is proportional to N, and which are well 
described by the BGS (see, for example, S. A. Cannas, A. C. N. Magalhaes, and F. 
A. Tamarit, Phys. Rev. B 61(2000) 11521). 

[5] Lists of many applications of the nonextensive statistics are available at URL: 
http://tsallis.cat.cbpf.br/biblio.htm. 

[6] T. Padmanabhan, Phys. Rep. 188 (1990) 285. 

[7] M. Schmidt, Phys. Rev. Lett. 86 (2001) 1191: see also, C. Tsallis, B. J. C. Cabral, 
A. Rapisarda, and V. Latora, e-print: cond-mat/0112266. 

[8] The specific heat in micro canonical ensembles of BGS may be not always positive 
(e.g. Ref. 6). 

[9] A. Lavagno, G. Kaniadakis, M. Rego-Monteiro, P. Quarati, and C. Tsallis, Astrophys. 
Lett. Commun. 35 (1998) 449. 

[10] A. Taruya and M. Sakagami, Physica A 322 (2003) 285. Phys. Rev. Lett. 90 (2003) 
181101: In the latter paper, it is shown that when Tsallis' old formalism is employed 
[Ref. 1], the q — n relation is given by n = l/(g — 1) + (3/2) leading to q = 1.05 ~ 1.66 
for n = 3 ~ 20, where n is the polytrope index. 

[11] M. L. Lyra and C. Tsallis, Phys. Rev. Lett. 80 (1998) 53. 

[12] G. Wilk and Z. Wlodarczyk, Phys. Rev. Lett. 84 (2000) 2770. 

[13] C. Beck, Phys. Rev. Lett. 87 (2001) 180601. 

[14] A. K. Rajagopal, C. S. Pande, and S. Abe, e-print: cond-mat/0403738. 

[15] S. Abe, S. Martinez, F. Pennini and A. Plastino, Phys. Lett. A 281 (2001) 126. 

19 



[16] S. Abe, Phys. Lett. A 263 (1999) 424: ibid. 267 (2000) 456 (erratum). 

[17] M. Sakagami and A. Taruya, Continuum Mechanics and Thermodynamics 16 (2004) 
279. 

[18] Y. Kakehashi, Adv. Phys. 53 (2004) 497; related references therein. 
[19] Y. Suezaki, Phys. Lett. 38A (1972) 293. 

[20] H. Shiba and P. A. Pincus, Phys. Rev. B 5 (1972) 1966: Prog. Theor. Phys. 48 
(1972) 2171. 

[21] U. Bernstein and P. Pincus, Phys. Rev. B 10 (1974) 3626. 

[22] M. E. Kozlov, V. A. Ivanov and K. Yakushi, J. Phys. C: Cond. Mat. 8 (1996) 1011. 
[23] H. Hasegawa, e-print: cond-mat/0408669. 

[24] M. Luban, J. Magn. Magn. Mat. 272-276 (2004) e635; H. Kachkachi and D. A. 
Garanin, e-print: cond/mat/0310694. 

[25] W. A. de Heer, P. Milani, and A. Chatelain, Phys. Rev. Lett. 65 (1990) 488. 

[26] J. P. Bucher, D. C. Douglass, and L. A. Bloomfield, Phys. Rev. Lett. 66 (1991) 3052. 

[27] S. E. Aspel, J. W. Emmert, J. Deng, and L. A. Bloomfield, Phys. Rev. Lett. 76 
(1996) 1441. 

[28] P. Gambardella, A. Dallmeyer, M. Malti, M. C. Malagoll, W. Eberhardt, K. Kern, 
and C. Carbone, Nature 416 (2002) 301. 

[29] Y. Yamamoto, T. Miura, H. Suzuki, N. Kawamura, T. Nakamura,K. Kobayashi, T. 
Teranishi, and H. Hori, Phys. Rev. Lett. 93, 116801 (2002). 

[30] T. Shinohara and T. Sato, Phys. Rev. Lett. 91 (2003) 197201. 

[31] A. V. Postnikov, M. Briiger, and J. Schnack, e-print: cond-mat/0404343. 

[32] D. Gatteshi, R. Sessoli, and A. Cornia, Chem. Commun. 9 (2000) 725. 

[33] H. Hasegawa, cond-mat/0501126. 



20 



[34] If the half-filled electron occupancy is satisfied for both de and d'-f subbands, the 
number of eigenvalues for the CE of the iV-site Hubbard model with the orbital 
degeneracy becomes C]$£ € C^^' y where D e = 2 and D 1 = 3 for 3d transition metals. 

[35] M. S. Reis, J. P. Araiijo, V. S. Amaral, E. K. Lenzi, and I. S. Oliveira, e-print, 
cond-mat/0207245. 

[36] S. Martinez, F. Pennini, and A. Plastino, Physica A 282 (2000) 193. 

[37] M. Portesi, A. Plastino, and C. Tsallis, Phys. Rev. E 52 (1995) R3317. 



21 



Figure 1: The temperature dependences of the energy E q of grand-canonical ensembles 
for h = with (a) U/t = 0, (b) 5 and (c) 10 calculated by the method A (GCE-A), and 
those with (d) U/t = 0, (e) 5 and (f) 10 calculated by the method B (GCE-B): q = 1.0 
(bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) and 2.0 
(solid curves). 

Figure 2: The temperature dependences of the entropy S q of grand-canonical ensembles 
for h = with (a) U/t = 0, (b) 5 and (c) 10 calculated by the method A (GCE-A) for 
h = 0, and those with (d) U/t = 0, (e) 5 and (f) 10 calculated by the method B (GCE-B): 
q = 1.0 (bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) 
and 2.0 (solid curves). 

Figure 3: The temperature dependences of the specific heat C q of grand-canonical ensem- 
bles for h = with (a) U/t = 0, (b) 5 and (c) 10 calculated by the method A (GCE-A), 
and those with (d) U/t = 0, (e) 5 and (f) 10 calculated by the method B (GCE-B): q = 1.0 
(bold solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) and 2.0 
(solid curves). 

Figure 4: The temperature dependences of the susceptibility Xq °f grand-canonical en- 
sembles for (a) U/t — 0, (b) 5 and (c) 10 calculated by the method A (GCE-A), and those 
for (d) U/t = 0, (e) 5 and (f) 10 calculated by the method B (GCE-B): q = 1.0 (bold 
solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) and 2.0 (solid 
curves) . 

Figure 5: The temperature dependence of 1/xq m the atomic limit (t/U — 0) of grand- 
canonical ensembles for various q values calculated by (a) the method A (GCE-A) and 
(b) the method B (GCE-B). The temperature dependence of T q (T) in the atomic limit 
(t/U = 0) for various q values calculated by (c) the method A and (d) B. The q dependence 
of the Curie constant T q (T) at T/U = and 10 calculated by (e) the method A and (f) 
B. Chain curves in (e) and (f) express r ? (oo) = g4 2 ^~ 1 ^ and r ? (oo) = q, respectively (see 
text). 

Figure 6: The temperature dependences of the specific heat C q of canonical ensembles for 
h — with (a) U/t = 0, (b) 5 and (c) 10 calculated by the method A (CE-A), and those 
with (d) U/t = 0, (e) 5 and (f) 10 calculated by the method B (CE-B): q = 1.0 (bold 
solid curves), 1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) and 2.0 (solid 
curves) . 

Figure 7: The temperature dependences of the susceptibility Xq °f canonical ensembles 
for (a) U/t = 0, (b) 5 and (c) 10 calculated by the method A (CE-A), and those for (d) 
U/t = 0, (e) 5 and (f) 10 calculated by the method B (CE-B): q = 1.0 (bold solid curves), 
1.1 (dotted curves), 1.2 (dashed curves), 1.5 (chain curves) and 2.0 (solid curves). 
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